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Abstract. In recent single-particle tracking (SPT) measurements on Listeria monocytogenes motil- 
ity in vitro, the actin-based stochastic dynamics of the bacterium movement is analyzed statistically 
(Kuo and McGrath, 2000). The mean-square displacement (MSD) of the detrended trajectory ex- 
hibit a linear behavior; it has been suggested that a corresponding analysis for the Brownian ratchet 
model (Peskin, Odell, & Oster, 1993) leads to a non-monotonic MSD. A simplified version of the 
Brownian ratchet, when its motion is limited by the bacterium movement, is proposed and analyzed 
stochastically. Analytical results for the simple model are obtained and statistical data analysis is 
investigated. The MSD of the stochastic bacterium movement is a quadratic function while the 
MSD for the detrended trajectory is shown to be linear. The mean velocity and effective diffu- 
sion constant of the propelled bacterium in the long-time limit, and the short-time relaxation are 
obtained from the MSD analysis. The MSD of the gap between actin and the bacterium exhibits 
an oscillatory behavior when there is a large resistant force from the bacterium. The stochastic 
model for actin-based motility is also mathematically equivalent to a model for integrate-and-firing 
neurons. Hence our mathematical results have applications in other biological problems. For com- 
parison, a continuous formalism of the BR model with great analytical simplicity is also studied. 
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1. Introduction 

Actin polymerization plays an important role in nonmuscle cell mechanics, motility, and func- 
tions (Pollard et al, 2000; Pantaloni et al., 2000). In recent years, quantitative analyses of the 
molecular mechanism for actin-based motility are made possible by both laboratory experiments 
on Listeria monocytogenes (see van Oudenaarden and Theriot, 1999, and the references cited within) 
and a series of insightful mathematical models (Hill, 1981, 1987, Peskin et al., 1993, Mogilner and 
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Oster 1996). The interaction between experimental observations and theoretical ideas has generated 
exciting research in biophysics and mathematical biology. 

Following the seminal work of Peskin et al. (1993), a sizable literature now exists on math- 
ematical models and analyses of the polymerization-based motility, known as Brownian ratchet 
(BR). Even though the original model on fluctuations is clearly a probabilistic one, it was cast 
mathematically in terms of the difference and differential equations with only a minimal stochas- 
tic interpretation. In the subsequent development, this stochastic nature of the model often has 
been obscured. In experimental laboratories, on the other hand, researchers often use Monte Carlo 
simulations to model the biological problem, partly because the data are inevitably stochastic. 

This situation has prevented a truly quantitative understanding of the actin-based motility 
and a closer interaction between the experimental measurements and mathematical modeling. In 
a recent experiment, Kuo and McGrath (2000) used the highly sensitive single-particle tracking 
(SPT) methodology to measure the stochastic movement of L. monocytogenes propelled by actin 
polymerization. The seemingly random data are then analyzed statistically in terms of the mean- 
square displacement (MSD). The exquisite data with nanometre precision reveals the discrete steps 
in the bacteria movement, presumably due to the actin polymerization, one G-actin monomer at a 
time. 

The stochastic nature of the BR, and the statistical treatment employed in experimental data 
analyses, necessitate a mathematical analysis of the BR model in fully stochastic terms. This is 
the main objective of the present work. Furthermore, Kuo and McGrath (2000) suggested that 
the BR movement, after detrending, exhibits a non-monotonic MSD. We shall investigate these 
practical issues as well. The significance of the stochastic interpretation is that one needs only to 
think about a single BR, and can derive theoretical MSD to compare with experiments. 

In order to clearly present the stochastic approach to the BR, we study only a special, but 
relevant, case of the generic BR model proposed by Peskin et al (1993). This restriction makes 
the model easily analyzed analytically. Interestingly, the mathematical model is also identical to 
one for integrate-and-firing neuron proposed many years ago by Gerstein and Mandelbrot (1964). 
In recent years, integrate-and-firing model has become one of the essential components in neural 
modeling (Hopfield and Herz, 1995). The stochastic model is quite basic; therefore we expect that 
our mathematical results also have applications in other branches of mathematical biology. The 
fractal nature of such model has also been discussed recently (Qian et al., 1999). 

All the mathematical background on stochastic processes used in this work can be found in the 
excellent text by Taylor and Karlin (1998). To help the readers who are not familiar with some 
of the stochastic mathematics, italic font is used for the key words when they first appear in the 
paper. 

2. Stochastic Formulation of a Brownian Ratchet Model 

(i) We consider an F-actin polymerizes in a 1-dimensional fashion with the rate of monomer 
addition a and the rate of depolymerization (3. a is a pseudo-first order rate constant which is 
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proportional to the G-actin monomer concentration. Each G-actin monomer has a size of 6. Hence 
the actin polymerization is modeled as a continuous-time random walk (Hill, 1987). We shall take 
the growing direction as positive, and denote the position of the tip of the actin filament by X(t) 
which is a stochastic process taking discrete values k6, where k is an integer. 

(ii) We assume that a bacterium is, in the front of the growing actin filament, located at Y(i): 
X(t) < Y(t). The bacterium has an intrinsic diffusion constant D^, and experiencing (or exerting) 
a resistant force F in the direction against the actin polymerization. In the absence of the actin 
filament, the bacterium movement is a Brownian motion with a constant drift rate —F/iji,. Since 
a bacterium is a living organism, the Db and the ijh are not necessarily related by the Einstein 
relation r]f,Di, = ksT for inert equilibrium objects. 

(iii) The F-actin and the bacterium interact only when they encounter: X(t) = Y(t). The actin 
filament, however, can not penetrate the bacteria wall. Therefore, the motion of the bacterium and 
the actin polymerization are coupled via a reflecting boundary condition at X(t) = Y(t). 

(i)-(iii) are the basic assumptions of the generic BR model first proposed by Peskin et al. (1993). 
In the present work, we shall further assume that (iv) the a is sufficiently large and (v) /? ~ 0. 
Therefore, whenever the gap A(t) = Y(t) — X(t) = 5, the gap will be immediately filled by a 
G-actin monomer, and the polymer does not depolymerize. These two assumptions correspond to 
a rapid polymerization condition under which the bacteria movement is the rate-limiting process 
in the overall kinetics. 

Fig. 1 shows the basic, stochastic behavior of X(t), Y(t), and A(t). Kuo and McGrath (2000) 
also introduced a detrended Y(t). Let v be the mean velocity of the bacterium movement Y(t), 
then the detrend Y(t) is defined as Y(t) = Y(t) — vt. 

Let be the time for incorporating the kt/i G-actin monomer. Then at time ^fc, X(^^) = 
Y(^fc) — When t > Y(t) follows a Brownian motion with diffusion constant D, drift rate 
—F/rjf), and reflecting boundary at kS. Y(t) moves stochastically and when it reaches {k + 1)6, 
denoted the time by ^k+i^ the (k-|-l)t/i G-actin monomer is incorporated. Then the process repeats. 
The waiting time for the next monomer to be incorporated is a random variable, we shall denote 
it by T: ^^^i = ^f. + T. This is our stochastic formalism for the BR model. Our analysis focuses 
on the stochastic properties of the random variable T. 

Fig. 2 shows the mean-square displacement (MSD) of the stochastic data in Fig. 1. The MSD 
for a stochastic processes X(t) is defined as 

MSD{t) = E[{X{r + t)-X{t))^] (1) 

which is a powerful analytical tool for analyzing stochastic processes with independent increments 
or stationarity. The E[. . .] in Eq. |l| denotes the expectation of random variables. For a stochastic 
process with independent increments, MSD(t) is further simplified into E [(X(t) — X(0))^] . In the 
case of a stationary process, its MSD is directly related to the correlation function: 

£;[X(t)X(0)] = E[X.^] - \mSD{t). (2) 
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The significance of MSD is that it can be obtained through a statistical analysis of stochastic 
experimental data (Qian et al., 1991). It is the essential link between the experimental mea- 
surements on fluctuations and stochastic mathematical models. For an experimental time series 
{xn\0 < n < N}, the MSD is defined as: 

^ N—m 
fc=0 

The statistical relation between the experimentally determined MSD in Eq. ^ and the theoretical 
MSD in Eq. || can be found in the paper by Qian et al. (1991). 

3. Basic Properties of the Model: Analytical Results 

Mean Waiting Time and Waiting Time Distribution. The time interval T between 
the repeated incorporation of successive actin monomer is the exit time of a diffusion process. By 
exit time T^, we mean the time a Brownian particle takes to reach 5 the first time, starting at z 
(0 < z < 5). Clearly is a random variable; its expectation T{z) = ^'[T^], known as mean first 
passage time, is the solution to the differential equation (Taylor and Karlin, 1998) 

D,T':, - (F/r?5)r; = -1 (4) 

with boundary conditions T^(0) = and T{5) = 0. Hence 

T{z) = !^ {e^^h.D, _ ,F./,,D,^ ^ ni.{z-5) _ 

Therefore, 

where uj = F5/{r]hDb) is the nondimensionalized resistant force. 

The probability density function /x^ (t) for the waiting time, can be obtained in terms of its 
Laplace transform, also known as the characteristic function of the random variable T^, QTiz,i^) 
~ /o°° /t^ (i)e~'^*dt which satisfies the following differential equation (Weiss, 1966) 

HI H = '^<9t(z, i^) (7) 

with boundary condition (?(5t(0, v)/dz = and Qt:{5, v) = \. Note Eq. ^ is a special case of Eq. 
I for T{z) = -dQT{z,0)/du. 

Eq. ^ can be analytically solved: 

where 
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Therefore, the variance in the waiting time 



5^ \ Se^^^ - (lOw - 6)e^ + - 2a; - 9 
Dl 



Var[T] = ( ^ ) -v---^- (9) 



If there is no resistant force from the bacteria, a; = and we have a simple expression 

Qt(0,i^) = (cosh^^V^) ^ (10) 

Renewal Processes, The Statistical Properties ofX.{t) and Y{t). With T as the waiting 
time, the tip of the rapid growing actin filament, X(t), is a renewal process. There is a large 
literature on this subject. The most relevant result to our model is the elementary renewal theorem 
for large t 

Em)] ^ (11) 

Therefore as a renewal process, a BR executes successive steps with size 6 and average time -E'[T]. 
The mean velocity of the BR, thus, is 



(12) 



E[T] \ 6 ) e^ -\-uj' 

This result is in agreement with that of Peskin et al. (1993). 

Furthermore from the theory of renewal process (Taylor and Karlin, 1998) 

Yar\m)\ « ^^l^t = o\ (13) 



where, according to Eqs. |6| and ^, 

i^*^ - (lOu; - 6 

\-uj) 



2 Se^'^ - (lOu; - 6)e'^ +a;2 - 2lj - 9 o 



The MSD for X(t), therefore, is 

E [(X(t) - X(0))2] pa ah + {vtf (15) 

which is a quadratic function of t. The expression for cr^ is a new result of the present work, which 
is comparable with experimental data. Fig. 3 shows the dependence of v and o"^ as functions of 
ijj = F5/r]i,Di,, the resistant force from the bacterium. If w = 0, then v = 2Di,/5 and o"^ = 14L';,/3. 

Since Y(t) — X(t) < 6 while both increase linearly with t, for large t Y(t) ~ X(t) with an 
error less than 6, the size of a single actin monomer. Strictly speaking, the Y(t) is not a stochastic 
process with independent increments. However, the error involved, again, is only on the order of 
the size of a single G-actin. To understand the statistical correlation of Y(t) within each "step", 
see the section below on the gap. 
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Detrend Y(f) and Its Statistical Properties. In the recent experimental work (Kuo and 
McGrath, 2000), the detrended Y(t) has also been reported, which can be defined as Y{t) = 
Y(t) — vt, where v is the mean velocity. For n5 < Y(t) < (n + 1)5, 

Y{t) ^ Y{t) -vt = Y{t) - Y{U +nd- v^^ - v{t - ^ = Y(t) - vt + n6 - 

(16) 

in which random variable ^„ is the time for Y{t) to reach n5 the first time, r = t — and v 



is given in Eq. [1^. Hence, < r < — = = T, with its expectation, variance, and 

characteristic function given in Eqs. |6|, |9|, and ^, respectively. 
The statistical properties of Y{t) are readily to be calculated: 



E 



Y{t) =E(Y{t)]-vt^O, (17) 



Var 



Y{t) = Var [Y(t)] + nvVar [T] ^ a'^t. (18) 



Thus, we see that the detrend Y{t) does not become stationary with increasing time. While its 
expectation is zero, its variance increases linear with the time t, the epitome of a symmetric random 



movement. The parameter ^ is the effective diffusion constant of the BR. 

Statistical Properties of the Gap. The gap between the bacteria, Y{t), and the tip of the 
actin filament X(t), A(t) = Y{t) — X(t) behaves completely different from the detrend Y(t). It 
reaches asymptotically to stationarity. 

We can provide a reasonable estimation for the relaxation time for the gap to reach its station- 
arity from the largest nonzero eigenvalues {/j.) of the diffusion operator 

Db:r^ + -—)u{x)=fiu{x) (19) 
dx^ rjh dx J 

under the boundary condition Dhu'{0) + {F/r]h)u{0) = u{6) = 0. See Appendix for details. All 
the eigenvalues are real and < 0, ^{z) = —^{z"^ +0;^), where the z are the roots of the tran- 

2 2 

scendental equation cos 2: = ^^^f^- Fig. 4 suggests that the largest eigenvalue corresponds to the 
exit time which increases with the resistant force. For resistant force F ^ Irj^Dh/S, there is a 
separation between the time scale for the exit and the time scale for establishing a quasi-stationary 
distribution for A(t) (Appendix). Fig. 5 shows the MSD for A(f), which is directly related to the 
correlation function for the stationary process (Eq. ^). After normalized by 2yar[A], the MSD 
are approximately the same for < u; < 6. The correlation time decreases with uj for uj = 2, 4, 6, 
and 12 (i.e., p = 0.55, 0.6, 0.65, and 0.8), corresponding to the second largest eigenvalue in Fig. 4, 
When there is a large resistant force F, the exit time T has a small relative variance (Eqs. |^ 
and ^ and the exit becomes an event with sufficient regularity. This is reflected in the oscillation 
of the MSD in Fig. 5. 

4. An Analytical Analysis of a Continuous Stochastic Formalism of BR 
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We can replace the assumptions (iv) and (v) in the Section 2 with a continuous model for 
the discrete polymerization. In other word, we approximate the random walk by a diffusion with 
diffusion constant and drift rate (Feller, 1957; Hill, 1987): 

Da = ia + (3)5^/2, Va = ia- 13)5. (20) 

where (3 and a are first-order and pseudo-first-order rate constants for the depolymerization and 
polymerization, 5 is the size of a G-actin monomer. The dynamic equation governing the probability 
density function (pdf) P-x^{x.,t) for the stochastic processes X(t) is: 

dt " dx'^ " dx 

where P-x{x,t) has the probabilistic meaning of Px{x)dx = Prob{x < X < x + dx}. Similarly, the 
dynamical equation for Brownian motion of the bacterium Y with resisting force F is, as before, 

dPY{y,t) ^^^ d^PYiy,t) ^ F dPYiy,t) 
dt 9y2 Qy 



These two equations are coupled since X < Y. We call Eqs K^, and 22 the continuous formalism 



for the BR. It represents a two-dimensional diffusion in the triangle region of x < y: 

dPKY{x,y,t) _ d^P^Y{x,y,t) d^P^Y{x,y,t) _ dP^Yix,y,t) F dPKY{x,y,t) 

ot ox'^ oy^ OX rjh oy ^^g-j 

The advantage of this version of the BR is its analytical simplicity. A coordinate transformation 
can be introduced: 

A = Y-X, Z=^°^ + "'^ (24) 

where A represents the gap between the tip of the actin filament and the bacterium, Z repre- 
sents an averaged position of X and Y, we shall call it the center of mass of the BR. With this 
transformation, the two differential equations are decoupled: 



dt \ a . 0, • y a < 

dPz{z,t) DgP^ f d^Pz{z,t) \ DbVa-DgF/r] f dPz{z,t) \ ^ ^ 

= TT-^ (— oo < z < +oo). 

dt Da + Dh\ dz^ J Da + Dh \ dz J ^ ' ^26) 

It can be immediately concluded from these two equations that the gap A(t) approaches to its 
stationary, exponential distribution (see Appendix) 

Z however increases steadily with an effective diffusion constant D^, = — "D^^' ^ mean 
velocity V, = ■ This result can be understood in terms of Newtonian mechanics: the 
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driving force from actin polymerization is Fa = tjaVa = ksTVa/Da, and the resistant force is F, 
and hence the net force on the BR is Fz = Fa — F with the frictional coefficient the BR (center of 
mass) being ksT/Dz- Hence 

' ksT Da + Db \Da kBT) Da + Db ' ^ ' 

The parameter Da and Va are defined in terms of the a, /3, and 5 in Eq. a is a pseudo-first order 
rate constant which is proportional to the G-actin concentration cq, as well as the probability of the 
gap A being greater than 5. Therefore, a is a function of external force F; it can be determined 
in a self-consistent manner by the transcendental equation: 



a{F) = ttQCQ / P^{s)ds = aocoexp 
Js 



{a-p)6 + F/rjb ' 
' {a + 13)5^/2 + Db . 



(29) 



where ag is the intrinsic, second-order rate constant for polymerization. We see that in Eq. 28 
is a linear function of resistant force F explicitly; however, nonlinearity arises since Da and Va are 
implicit functions of the resistant force F, via a{F). In other words, the resistant force F slows 
down the BR by two different mechanisms: a linear Newtonian resistance and also a reduction in 
the rate of polymerization via a diminished gap. Eq. ^ has the general form Oi{F) oc g,~'''P^/^BT 
with an entropic barrier. 

Vz in Eq. ^is necessarily smaller than Va, indicating that the bacterium retards the polymeriza- 

( nbDb\ 2(«-/3) 



tion. When F = rjbDbVa/Da = y^^^j a+fd ' bacterium completely stalls the polymerization. 

This yields the critical stalling force which agrees with the well known result of Hill (1987) for 

an inert object: [ksT / 5)ln[a/ (3). Furthermore, if we note that rjbDb = ksT, then the rate of 

polymerization against a resisting force F is at its maximal Vz = Va — FDa/ksT when Db oo. 

This is a result of Peskin et al. (1993) who first elucidated the crucial role played by the fluctuating 

"barrier" in the filamental growth. In a more general context, the dynamic characteristics of the 

"force transducer" by which the resisting force is applied to the growing tip of the filament is an 

integral part of the molecular process .[| 

We now show that in the limit of 5 ^ 0, our Vz from the continuous model is equivalent to the 

ratchet velocity derived by Peskin, Odell, and Oster (POO, 1993). Note that in the previous work 

the definition for the ratchet velocity was rather convoluted; The Vz in the present model is more 

'^In the work of Hill (1981), this issue is not considered because of its quasi-thermodynamic approach. The 
applied force was assumed to have an instantaneous dynamic characteristics. From equilibrium thermodynamics 
one knows that polymerization under resisting force F has a{F) / l3{F) = e~^^^''^^ . So how does F contribute 
to a{F) and P{F) individually? A splitting parameter r is defined as a{F) — a{0)e~'^^^''''^^ and consequently 
I3{F) = /3(0)e('-'"'-^*/'=«^, where a(0) and /?(0) are the a and P in the main text. The rate of polymerization under 
force F, therefore, is 

[a{F) - I3{F)]S = (^Q(0)e""^*'"''«^ - /3(0)e''""'-^*'"=«^) 5. (30) 

With very smaU 5, this gives [q(0) - /3(0)]5 - FS'^[r q(0) + (1 - r)/3(0)]/fcsT which should be compared with K 
= (a ~ P)5 - F5^{a + l3)/2kBT from the main text. This indicates that the BR has a splitting factor of r = 1/2. 
However, Hill's analysis does not have the contribution from the dynamic characteristics of the barrier, i.e., Db- Eq. 
^ is the starting point of the recent work of Kolomeisky and Fisher (2001). 
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straightforward. From Eq. 2C we have a. 
these two expressions into 



we have 



hm Vr, 

5^0 



poo 



poo 



a 



^) and /? = i ( 

roc roo 

/ P^{A)dA-l3 P^iA)dA 
Js Jo 



-^ + ^). Substituting 



P^{A)dA = Va-DaPA{0) 



Da + Dt, 



In the derivation we have used Eq. 27. Note that the basic molecular parameters for polymerization, 
a, (3 and 5 are actually contained in the parameter Da and Va- In the mathematical limit of 5 — > 0, 
there is at the same time a and /? — > oo such that Da and Va are finite (Feller, 1957). 



5. Discussion 

Nanometre precision measurements on L. monocytogenes movement (Kuo and McGarth, 2000) 
have shown that the bacteria move with steps. Considering there are many actin filaments in a 
bundle which propels a bacterium, this observation indicates a synchronized filamental growth in 
the bundle. The synchronization is not inconsistent with a bundle of actin filaments propelling a 
bacterium with sufficiently small Brownian movement (i.e., small Dh). The significantly reduced 
Brownian motion is indeed observed experimentally, both in the direction parallel and perpendicular 
to the actin growth. 

There could be several explanations for the small Di,. a) Kuo and McGrath (2000) suggested 
an association between the bacterium and the actin structure, which leads to endorsing the two- 
dimensional BR with bending (Mogilner and Oster, 1996). b) It should be noted, however, that 
association-dissociation can also be introduced into the one-dimensional BR in the form of an 
attractive force between the actin and the bacterium; thus a nonzero F as function of {y — x). This, 
we suspect, will also lead to a reduced apparent D\, on a longer time scale, c) As we have pointed out, 
the Db of a living bacterium is not necessarily related to its physical size and frictional coefficient r/f,. 
A bacterium could have an internal mechanism, by utilizing its biochemical free energy, to localize 
itself near the tip of the actin filament with diminished Brownian motion. Finally, all existing 
models on BR have only dealt with single filaments. The continuous formalism we proposed here 
is in fact our initial step to extend the BR to a filamentous bundle. All these topics are currently 
under investigation. 

In a special Science issue on Movement: Molecular to Robotic, two articles reviewed recent 
progress on force and motion generated on the molecular level by two completely different bio- 
logical systems: motor protein movement and cytoskeletal filamental polymerization (Vale and 
Milligan, 2000; Mahadevan and Matsudaira, 2000). Both systems can move against resistant force 
by utilizing chemical free energy. In the abstract of the second article, it was stated "Not all biolog- 
ical movements are caused by molecular motors sliding along filaments or tubules. Just as springs 
and ratchets can store or release energy and rectify motion in physical systems, their analogs can 
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perform similar functions in biological systems." While there has been much work done on motor 
proteins and protein polymerization in connection to various cellular phenomena such as motility, 
less has been discussed about the fundamental physical principles of these two processes. It turns 
out that both molecular processes have a single, unified mathematical model which accounts for 
their chemomechanical energy transduction. 

Theoretical formalism for motor proteins are now well established (sec Jiilicher et al., 1997; 
Qian, 2000b, and references cited within). Since the motion of a single motor protein is Brownian, 
it has to be characterized in terms of probability distribution. The simplest model is that of Huxley 
(1957). This model corresponds to one on polymerization with nucleotide hydrolysis proposed by 
Dogtcrom and Lciblcr (1993). Both models addressed the important issue of nucleotide hydrolysis, 
but neglected the stochastic nature in the movement of motor protein and actin polymerization, 
respectively. Without the diffusion term, such mathematical model is known as random evolution 
(Pinsky, 1991). 

With the ATP cap and hydrolysis, the model for the stochastic dynamics of actin polymer- 
ization will be precisely in the same class of the models for single motor proteins (S.-D. Liang, 
G. Martinez, G.M. Odell, and H. Qian, work in progress). The experimental measurements on 
both systems also proceed with parallel paths, as demonstrated by Dogterom and Yurke (1997), 
and more recently Kuo and McGrath (2000). Similar to the measurements on load- velocity curves 
for motor proteins, Dogterom and Yurke measured the velocity as a function of resistant force for 
single microtubules growing in vitro. Analysis of their data suggests that under the stalled (crit- 
ical) condition, polymerization is in a noncquilibrium steady-state rather than a thermodynamic 
equilibrium (Kolomcisky and Fisher, 2001; Hill, 1987). 

All these experimental evidences indicate that the class of BR model (or augmented Huxley 
model) is a fundamental mathematical model for chemomechanical energy transduction. Recent 
work on the noncquilibrium statistical mechanics and thermodynamics of BR, in the context for 
single macromolecules in aqueous solution, also provided the mathematical model with a solid 
foundation in statistical physics of Boltzmann, Gibbs, and Onsager (Qian, 1998, 2000b, 2001a,b,c). 
The biological systems discussed in the two Science articles (Vale and Milligan, 2000; Mahadevan 
and Matsudaira, 2000) and the mechanistic, molecular models proposed are completely different. 
Yet they share fundamentally the same physiochemical principle which units both models in a 
quantitative fashion. The mathematical model seems to capture the basic principle for molecular 
movements and forces in cell biology. 

To mathematical biologists, BR is a class of models which is based on a similar physical model 
but can have many different mathematical representations and different degree of approximations. 
We have shown two such analyses in the present work. The essential feature of all the models can, 
and should, be presented in terms of their MSD, which provides the BR, in steady-state, with an 
effective diffusion constant(D^) and a mean velocity (Vz), both as functions of the resistant force. 
More subtle differences between models can be found in the transient behavior. The comparison 
between our analyses is summarized in the Table, in which the effective diffusion constant for the 



10 



continuous model 



(a + /3)(52 A/2 



A + (a + /3)(572 



when a — > oo, and the BR velocity 



Db{a - 13)6 - Dbja + p)Sio/2 



{Dh/S){2-oj) 



is a linear force-velocity relationship. 



discrete model 



3e^"-(10a)-6)e"+u;^-2a;-9 2 
(e'^-l-a;)3 




continuous model 



1 



2-LO 



It is seen that in the continuous model, the effective diffusion constant is always less than 
the Db, while in the discrete model, the = is a function of the resistant force. When the 

force is small, can in fact be greater than Db- This is a type of facilitated diffusion. 

It is important to point out that the results from our discrete analysis is invalid when the 
resistant F is sufficiently large, when the polymerization is near its stalling force. This is due to 
the assumption of infinite large a. This explains why there is no critical force in Fig. 3, at which 
the velocity v = 0. The more realistic model with finite a and /? does lead to a finite, positive 
stalling force (Peskin et al., 1993; Kolomeisky and Fisher, 2001). The valid regime for our discrete 
model is a rapid growing actin filament with the bacteria viscous drag being the limiting factor in 
the overall BR movement. 

Finally, it is worth pointing out that mechanical studies of cellular properties and functions 
can be approximately classified as for passive and active materials. The former can be understood 
in terms of the theories of viscoelasticity and polymer dynamics, see Qian (2000a) for a general 
approach to the problem. Materials with chemomechanical energy transduction are active. The 
fundamental difference is the nucleotide hydrolysis which leads to an irreversible thermodynamic 
nonequilibrium steady-state, with heat dissipation (Qian, 2001c), in the latter rather than the 
usual equilibrium. Thus, the BR model is also a natural generalization of the standard polymer 
theory for passive materials (Doi and Edward, 1986) to active materials for which T.L. Hill (1987) 
has coined the term "steady-state polymer" . There is a continuous intellectual thread in all these 
mathematical theories. 
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8. Appendix 



Diffusion with Drift in Semi-infinite Space with Noflux Boundary. To understand 
the dynamics of the gap in the continuous model, one needs to solve the time-dependent diffusion 



equation, Eq. |25|. In nondimensionalized form: 

ut = Uxx+i^Ux, (2; > 0), 



(31) 



with boundary condition Ux + uju = at x = and x = 00. Amazingly, classic texts on diffusion 
(Carslaw and Jaeger, 1959; Crank, 1975) did not give an explicit solution to the problem. Becasue 
of its central importance in the theory of BR, we give some explicit results below. 

The eigenfunction of the problem associated with the eigenvalue /j,{z) = —j{uj'^ + z^) is 

1 



u{x, t; z) 



z cos 



zx 



u! sm 



^)] e~^+^(^)*, (z>0). 



(32) 



and for ^ = 0, y/uje . Note there is a gap in ^ between = and the continuous spectrum 
/X < — The Sturm-Liouville eigenvalue problem has a complete orthonormal set 

ti(x, 0; z)n(x, 0, z)e''''dx = 5{z - z'). 



Therefore, the solution to Eq. ^ with initial data 5{x) is 

zdz 
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(33) 



which approaches to the exponential distribution uje when t ^ 00. From Eq. ^ we have 

n(x, t)dx = 1, 



xm(x, t)dx 
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e 



UJ t 
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dz 
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UJ 



7r(u;2 + 22)2 

(1 + 2T)er/ (^) - 2r + 2 



(34) 



where r = uP't/A. The curve in the square bracket is a universal curve, we shall denote it by 
1 — gap{T). See Fig. 6. It is interesting to point out that there is a sharp transition between an 
unlimited growth of x for lij < and a stationary state for x when uj > 0. This mathematical result 
is similar to that of Dogterom and Leibler (1993). Finally, 



/ x'^u{x,t)dx = / gap{s)ds 
Jo Jo 



The Dynamics of Gap A(t) in the Discrete Model. Eq. |19| has a second boundary 
condition u = at x = 6, which corresponds to nondimensionalized Eq. ^ with < x < 1 and 
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n(l) = 0. It has a set of discrete eigevalues. The eigenfunctions to the nondimensionahzed Eq. 
|l9| with eigenvalue fi{z) = —^{uj"^ + z^) < 0, are still given in Eq. but the z's are now the 

2 2 

discrete roots of the transcendental equation cos 2; = . When u) < 2, the equation for z has 

only real roots; hence the largest eigenvalues is /i < — If, however, co > 2, then there is a pair 
of imaginary ±iz*, \z*\ < co. Then the largest eigenvalue is = — (w^ — (z*)'^^ /4. Fig. 4 shows 
how the largest eigenvalue (smallest in magnitude) changes as functions of w. It is seen that the 
largest eigenvalue can be well represented by 1/£'[T] 



ui ^ or zi=uj\ ; 1. (35) 

The solution to the time-dependent Eq. |19| with initial data 5{x) can be obtained in terms of the 
u{x,t]zys in Eq. 32, and approximated by the first term: 



u{x,t) K, ——^ — 2T cos ( ) — w sin ( ) exp ( — — ^ t ] . (36) 



ZlX\ . ( Z\X 



2 ) V 2 



UJ uj'^ + zf 



2 



7r(w2 _^ ^2 

Therefore, 

' / X , 2zi ^ ^"^+"1 , . fzi 

u[xA)ax = — TT-e 2 4 gm — 

^ ^ 7r(cj2 + zl) V 2 

which tends to zero because the exit probability. However, the remaining probability for A(t) 
quickly approaches to a quasi-stationary distribution 



2sm(^j V 2 



(jp + zl . fZl 

sm — 

2zi V 2 



(l-x))exp ■ 



For large a;, this distribution approaches to the exponential distribution oje as expected (data 
not shown). 
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Figure 1: A set of examples, from Monte Carlo simulations, for the stochastic trajectories of Y(t), 
the movement of the bacterium, X(f), the movement of the tip of the actin filament, A(t) = 
Y(t) — X(t), the gap, and Y(t), the dctrended Y(t). Among the four types of data, only the 
A(i) approaches stationarity. Y(i) has zero expectation but with a linear MSD, a characteristic of 
Brownian motion without drift. Both X(t) and Y(t) show the typical diffusion with a drift (Qian 
et al., 1990). 
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Figure 2: The MSD calculated for the four types of data in Fig. 1. Y(t), the movement of the 
bacterium, X(t), the movement of the tip of the actin filament, A(t), the gap, and the detrend 
Y(t). As expected, after a brief period of time, the Y(t) and X(t) are almost indistinguishable; the 
gap between them quickly reaches stationarity. The detrend Y{t) shows a linear MSD, as observed 
in the experiments (Kuo and McGarth, 2000). 
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Figure 3: The velocity vS/D}, and MSD /Di^, nondimensionalized, as functions of the resistant 
force on bacterium, oj = F5/r]i,Di,. In the SPT experiments, the velocity can be obtained as the 
quadratic term in the MSD of Y(t), the can be obtained either as the linear term in the MSD 
of Y(t), or the MSD of the detrend Y(t). The velocity v is given in Eq. 12 and the a is given in 
Eq. 11. 
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Figure 4: Numerical computation for the three smallest eigenvalues (in magnitude, all eigenvalues 
are negative) of Eq. |l^, as function of the resistant force uj. These are the most relevant 

modes in the relaxation (and correlation function) of the gap, A(t), approaching to stationarity. 
The smallest eigenvalue corresponds to the exit time, shown in the figure (labeled 1/-E[T]). 
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Figure 5: The normalized MSD for the A(t), from Monte Carlo simulation, for different resis- 
tant force F which is related to the probability (p) shown in the figure: F/(t]i,5) = 0.1(2p — 1). 
A standard MSD for a stationary process is directly related to its time correlation function 
2 (£'[A^(t)] - E[A{t)A{0)]), with its asymptote being the 2Var[A] when t ^ oo. In the simula- 
tions, the diffusion constant is Db/S"^ = 0.005. < A{t) < 1; the stationary E[A] = 0.33, 0.27, 0.21, 
0.16, and 0.08 for p = 0.5 — 0.8 respectively; the corresponding relative variances Var[A\/E'^[A] = 
0.53, 0.64, 0.79, 0.91, and 0.96. The relative variance for an exponential distribution is 1. 
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Figure 6: The function gap{t), defined in Eq. 34, is the time course for the the mean gap size to 
approach to its stationarity. For a; > 0, the mean gap size approaches to a finite size, while for 
Lu < 0, the mean gap size grows without bound. 
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